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Spectroscopy of B c Mesons in the Relativized Quark Model 
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We calculate the spectrum of the charm-beauty mesons using the relativized quark model. Using 
the wavefunctions from this model we compute the radiative widths of excited cb states. The 
hadronic transition rates between cb states are estimated using the Kuang-Yan approach and are 
combined with the radiative widths to give estimates of the relative branching ratios. These results 
are combined with production rates at the Tevatron and the LHC to suggest promising signals 
for excited B c states. Our results are compared with other models to gauge the reliability of the 
predictions and point out differences. 
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I. INTRODUCTION 

The charm-beauty (B c ) quarkonium states provide a 
unique window into heavy quark dynamics and there- 
fore an important test of quantum chromodynamics. Al- 
though they are intermediate to the charmonium and 
bottomonium systems the properties of B c mesons are 
a special case in quarkonium spectroscopy as they are 
the only quarkonia consisting of heavy quarks with dif- 
ferent flavours. Because they carry flavour they cannot 
annihilate into gluons so are more stable with widths less 
than a hundred keV. Excited B c states lying below BD 
(and BD* or B*D) threshold can only undergo radia- 
tive or hadronic transitions to the ground state B c which 
then decays weakly. This results in a rich spectroscopy 
of narrow radial and orbital excitations below p(*)£)(*) 
threshold which are more stable than their charmonium 
and bottomonium analogues: There are two sets of S- 
wave states, as many as two P-wave multiplets (the IP 
and some or all of the 2P) and one D-wave multiplet be- 
low BD threshold. As well, the F-wave multiplet is suffi- 
ciently close to threshold that they may also be relatively 
narrow due to angular momentum barrier suppression of 
the Zweig allowed strong decays. 

The discovery of the B c meson by the Collider De- 
tector at Fermilab (CDF) Collaboration l] in pp col- 
lisions at y/s — 1.8 TeV has demonstrated the pos- 
sibility of the experimental study of this system and 
has stimulated considerable interest in B r spectroscopy 

asiiiaaiiiamEEiEiElm. c a i- 

culations of B c cross sections at hadron colliders predict 
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that large samples of B c states should be produced at 
the Tevatron and at the LHC opening u p this new spec- 
troscopy in h m m m m 11 hm m eu . At the 

Tevatron it is estimated that O(10 7 ) B c mesons should be 
produced for 1 fb _1 of integrated luminosity while at the 
LHC O(10 9 ) B c mesons are expected to be produced for 
L = 100 fb _1 . These numbers are highly sensitive to the 
Pt and rapidity cuts used to extract the signal [l!J |2(| • 
The B* cross sections are expected to be 50-100% larger 
than the B c cross sections Qja Wi, 123 and Cheung and 
Yuan [22, |25| predict excited P-waves to contribute 20% 
of inclusive B c production while D-wave states are ex- 

Pted to contribute about 2% J2g. Chang and Chen 
and Cheung [2l| estimate that the 2S states will be 
produced in roughly the ratio of 2S/1S ~ 0.6. It should 
therefore be possible to start exploring cb spectroscopy 
at the Tevatron, producing the IP and 25 states and 
possibly even some ID and 2P states with sufficient rate 
to be observed. At the LHC, with its higher luminos- 
ity, the D-wave cb states should be produce in a sizable 
number so that the LHC should allow the study of the 
spectroscopy and decay of B c mesons. 

In this paper we study the spectroscopy, including ra- 
diative transitions, of charm-beauty mesons using the rel- 
ativized quark model [&|23,|2ij. The model includes one- 
gluon exchange with a running coupling constant and a 
linear confining potential. It uses relativistic kinemat- 
ics and momentum dependent and nonlocal interactions. 
Although this model is not a rigorous calculation from 
first principles it gives a good account of most known 
mesons and baryons with only a few free parameters 
@ IH, |2Jj, |3(J, ill so that it provides a useful guide to 
missing states. We compare its predictions to those of 
other models with the aim of highlighting which predic- 
tions are most sensitive to details of the models and there- 



fore the most useful in distinguishing models. However, 
we are also interested in pointing out which predictions 
give the greatest agreement between models and there- 
fore offer the most robust signatures for experiments to 
look for. Observation of these states and measurement 
of their properties would provide valuable information 
distinguishing details of the various models. 

We start with a brief outline of the relativized quark 
model and comment on its similarities and differences 
with other quark model calculations. The spectroscopic 
predictions are given and compared to those of other 
models in the literature. This is followed by predic- 
tions for El and Ml radiative transitions and estimates of 
hadronic transitions based on the Kuang-Yan approach 
[32. l33j . We summarize existing predictions for some of 
the more prominent weak decays of the B c ground state 
as final states are an important ingredient in reconstruct- 
ing the B c mesons. In addition, leptonic decays measure 
the wavefunction at the ori gin and are therefore an ad- 
ditional test of the model J3J. We end by discussing 
some strategies for searching for excited B c mesons and 
studying their spectroscopy. 



II. SPECTROSCOPY 

In this section we give the mass predictions of the rel- 
ativized quark model [a, |2i| |23 for the charm-beauty 
mesons and compare those predictions with the predic- 
tions of other calculations. This model has ingredients 
common to many quark potential models [4l8ll9l [Kl[l2| . 
Almost all such models are based on some variant of the 
Coulomb plus linear potential expected from QCD. An 
interesting observation is that all recent models have ar- 
rived at the same slope for the linear confining potential 
of ~ 0.18 GeV 2 . Most models, as does ours, also in- 
clude the running constant of QCD, a s (Q 2 ). And finally, 
relativistic effects are often included at some level. The 
relativized quark model has been reasonably successful in 
describing most known mesons. Although cracks have re- 



cently appeared with the discovery of the D S j |35ll3til37| 
and A (3872) states [38J, these point to the need to in- 
clude physics which has hitherto been neglected such as 
coupled channel effects 39] . 

In the relativized quark model mesons are approxi- 
mated by the qq sector of Fock space, in effect integrat- 
ing out the degrees of freedom below some distance scale, 
/i _1 . This results in an effective potential, V(p, r), whose 
dynamics are governed by a Lorentz vector one-gluon- 
exchange interaction at short distance and a Lorentz 
scalar linear confining interaction. The basic equation 
of the model is the rest frame Schrodinger-type equation 



H\^) = [H + V q ^f)M=E\i>) 



where 



H = ■ /- 2 



(1) 



(2) 



The effective quark- antiquark potential, V qq (p, r), was 
found by equating the scattering amplitude of free 
quarks, using a scattering kernel with the desired Dirac 
structure, with the effects between bound quarks inside 
a hadron |4ll |42| . Due to relativistic effects the potential 
is momentum dependent in addition to being co-ordinate 
dependent. To first order in (v/c) 2 , V q q(p,r) reduces to 
the standard non-rclativistic result: 



V q -M r) - V(r) = H™ nf + Hr nt 
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where 



HH 6 r 



(3) 



(4) 



includes the spin-independent linear confinement and 
Coulomb-like interaction, 



Tjcont _ ^27r a s (r) ^ ^ j3/ 
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is the colour contact interaction, 
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is the colour tensor interaction, 

TTS.O. _ TT S -°-( Cm ) 

n qq — n qq 

is the spin-orbit interaction with 
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its colour magnetic piece arising from one-gluon exchange 
and 
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qq 
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the Thomas precession term. In these formulae a s {r) is 
the running coupling constant of QCD. 

To relativize the qq potential, the full Dirac scattering 
amplitude was used as a starting point which for on-shell 
qq scattering is exact. However for a strongly interacting 
system there will in general be off-shell behavior which 
we did not consider in addition to other simplifications 
such as neglecting more complex components of Fock 
space. We therefore built a semiquantitative model of 
relativistic effects by smearing the co-ordinate r over the 
distances of the order of the inverse quark mass by con- 
volving the potential with a Gaussian form factor and 
replacing factors of m~ with, roughly speaking, factors 
of (p 2 + m 2 )~ 1 > 2 . The details of this relativization proce- 
dure and the method of solution can be found in Ref. [y] . 
It should be kept in mind that because we neglected cou- 
pled channel effects and the crudeness of the relativiza- 
tion procedure we do not expect the mass predictions to 
be accurate to better than ~ 10 — 20 MeV. 



For the case of a quark and antiquark of unequal mass 
charge conjugation parity is no longer a good quantum 
number so that states with different total spins but with 
the same total angular momentum, such as the 3 Pi — 1 P\ 
and 3 Z?2 ~ l D2 pairs, can mix via the spin orbit inter- 
action or some other mechanism. Eqns. [8] and [!|] can be 
rewritten to explicitely give the antisymmetric spin-orbit 
mixing term: 



H-= + l(i?i- b -- 



1 



4 \3r 
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where 5_ = Sq — Sq. Consequently, the physical j = 1 
P-wave states are linear combinations of 3 Pi and l P\ 
which we describe by: 



P' =* Pi cos 9 nP + 3 Pi sin e nP 
P = - 1 Pi sin 8 nP + 3 Pi cos nP 



11 



with analogous notation for the corresponding L = D, 
P, etc. pairs. In eqn. ^] P = L = 1 designates the 
relative angular momentum of the QQ pair and the sub- 
script J — 1 is the total angular momentum of the QQ 
pair which is equal to L. Our notation implicitely im- 
plies L — S coupling between the quark spins and the 
relative orbital angular momentum. In the heavy quark 
limit in which the heavy quark mass toq — > 00, the states 
can be described by the total angular momentum of the 
light quark, j, which couples to the spin of the heavy 
quark and corresponds to j — j coupling. This limit gives 
rise to two doublets, one with j — 1/2 and the other 
with j = 3/2 and corresponds to two physically inde- 
pendent mixing angles 9 = — tan~ 1 (v / 2) — —54.7° and 
6 = tan -1 (1/a/2) ^ 35.3° HJ|. Some authors prefer to 
use the j — j basis |2J but since we solve our Hamiltonian 
equations assuming L — S eigenstates and then include 
the LS mixing we use the notation of eqn. ^] It is 
straightforward to transform between the L~S basis and 
the j — j basis. It will turn out that radiative transitions 
are particularly sensitive to the 3 L L - 1 L L mixing an- 
gle with predictions from different models in some cases 
giving radically different results. We also note that the 
definition of the mixing angles are frought with ambigui- 
ties. For example, charge conjugating cb into be flips the 
sign of the angle and the phase convention depends on 
the order of coupling L, Sq and Sq 43]. 

The Hamiltonian problem was solved using the follow- 
ing parameters: the slope of the linear confining poten- 
tial is 0.18 GeV 2 , m c = 1.628 GeV, and m b = 4.977 GeV. 
The predictions of our model are given in Fig. 1 and are 



compared to the predictions of other calculations in Ta- 
ble [I] Because the mixing angles defined in eqn. 1111 are 
important for predictions of radiative transitions those 
predictions are also given in Table [I] Although I have 
attempted to consistently give the masses and the mix- 
ing angles of the predicted eigenstates in the convention 
of eqn. 111! because not all authors have unambiguously 
defined their phase conventions I cannot guaranty that 
these results are free of inconsistantcies. 

The different models are in remarkable agreement with 
the differences, for the most part, within the expected 
accuracy of the models. This almost certainly indicates 
how the various models have converged to using simi- 
lar confining potentials and including a strong running 
coupling constant in the Coulomb piece of the poten- 
tial. The only significant difference is the larger spread 
(~ 70 MeV) for the ID multiplet centre of gravity pre- 
dictions. The spin-dependent splittings are also in rea- 
sonable agreement. Potential models can therefore be 
used as a reliable guide in searching for the B c excited 
states. An important difference in the predictions is that 
in the Eichten-Quigg calculation [2] the lPi states are al- 
most pure 3 Pi and l P\ with little mixing while in other 
models there is significant mixing. This arises from the 
much smaller expectation value of the off-diagonal mix- 
ing term (eqn. HU|) in the Eichten-Quigg calculation Q 
compared to the other models. Since, after rotating from 
the j — j basis to the L — S basis, the L — S mixing 
term given by Eichten and Quigg is in agreement with 
eqn. 1101 the differences in expectation values can only be 
attributed to differences in the cancellations between the 
short and long distance pieces in eqn. 1101 i.e. between 
the |^f and £ pieces. This reflects subtle differences in 
the qq potentials of different models. Because the El ra- 
diative transitions are sensitive to the 3 Pi — Pi mixing, 
the measurement of radiative transitions can be used to 
distinguish between the different models. The study of 
the B c spectroscopy will help test and refine the quark 
potential models but more importantly will test Lattice 
QCD, NRQCD, and pNRQCD etc. which are more di- 
rectly connected to QCD. 



III. RADIATIVE TRANSITIONS 

Radiative transitions will likely play an important role 
in the discovery and identification of B c states. In this 
section we calculate the El and Ml radiative widths. 
The partial width for an El radiative transition between 
states in the nonrelativistic quark model is given by [45j 



J 



, 2S+1 T , i2S' + 1t 1 



f2S' + lj I I „ |„2S+1 T v |2 



T(n T Lj -> n Lj, +7) = -(e Q ) au C fi d S s> \ (n T L 7 , | r \ n T L, 7 ; 



(12) 



r 



where 



m b e c - m c e- b 
rrn, + m c 



(13) 



e c = 2/3 is the c-quark charge and e b = —1/3 is the 
6-quark charge in units of |e|, a is the fine-structure con- 



TABLE I: Predicted masses and Spin-Orbit mixing angles. The first column labelled GI is the present work. The P[ 
D' 2 — D 2 , and P3 — F3 states and mixing angles are defined using the convention of eqn. 1111 



Pi, 



State 


GI[6] 


EFG [12] 


FUII [10] 


GKLT [5] EQ [2] 


GJ [91 


ZVR [8] 


Lattice [44]" 


l 3 Si 


6338 


6332 


6341 


6317 


6337 


6308 


6340 


6321 ± 20 


1^0 


6271 


6270 


6286 


6253 


6264 


6247 


6260 


6280 ± 30 ± 190 


1 3 P 2 


6768 


6762 


6772 


6743 


6747 


6773 


6760 


6783 ± 30 


1P{ 


6750 


6749 


6760 


6729 


6736 


6757 


6740 


6765 ± 30 


lPi 


6741 


6734 


6737 


6717 


6730 


6738 


6730 


6743 ± 30 


l 3 Po 


6706 


6699 


6701 


6683 


6700 


6689 


6680 


6727 ± 30 


e-LP 


22.4° 


20.4° 


28.5° 


17.1° 


~-2° 


25.6° 




33.4 ±1.5° 


2 3 Si 


6887 


6881 


6914 


6902 


6899 


6886 


6900 


6990 ± 80 


2 1 S 


6855 


6835 


6882 


6867 


6856 


6853 


6850 


6960 ± 80 


2 3 P 2 


7164 


7156 




7134 


7153 




7160 




2Pi 


7150 


7145 




7124 


7142 




7150 




2Pi 


7145 


7126 




7113 


7135 




7140 




2 3 Po 


7122 


7091 




7088 


7108 




7100 




@2P 


18.9° 


23.2° 




21.8° 


-17° 








3 3 S! 


7272 


7235 






7280 




7280 




3^0 


7250 


7193 






7244 




7240 




1 3 P>3 


7045 


7081 


7032 


7007 


7005 




7040 




1P>2 


7036 


7079 


7028 


7016 


7009 6 




7030 




1P>2 


7041 


7077 


7028 


7001 


7012 6 




7020 




l 3 Di 


7028 


7072 


7019 


7008 


7012 




7010 




OlD 


44.5° 


-35.9° 




-34.5° 










1 3 P 4 


7271 












7250 




1*3 


7266 












7250 




1P3 


7276 












7240 




1 3 F 2 


7269 












7240 




OlF 


41.4° 

















"The error estimates are taken from Ref. 111! . 
6 We identify the 1D 2 and 1D 2 states with the 1 D 2 and 3 D 2 states 
of Eichten and Quigg 2j. 



stant, u) is the photon's energy, and C/j is given by 



(14) 



For convenience the Cfi coefficients are listed in Table 

HH The matrix elements (n' 2S ' +1 L[,,\ r \n 2S+1 L.j) are 
given in Table ITT1 and were evaluated using the wavefunc- 
tions given by the relativized quark model [6J . Relativis- 
tic corrections are implicitly included in these El transi- 
tions through Siegert's theorem Ha, E3, El] , by including 
spin-dependent interactions in the Hamiltonian used to 
calculate the meson masses and wavef unctions. The El 
radiative widths are given in Table [n] and compared to 
other predictions in Table HVI 

Most of the predictions for El transitions are in quali- 
tative agreement. While most differences are due to dif- 
ferences in phase space arising from different mass pre- 
dictions the more interesting differences arise from wave- 
function effects. The largest differences are for decays 
involving the Pi and P[ states which are mixtures of the 



spin singlet l P\ and spin triplet 3 Pi states. These can be 
traced back to the different 3 Pi — 1 Pi mixing angles pre- 
dicted by the different models. Wavefunction effects also 
appear in decays from radially excited states to ground 
state mesons such as 2 3 Pq — > l 3 Si 7 which varies from 1 
to 22 keV. Because the 2 3 Po has a node in its wavefunc- 
tion there will be a cancellation between different pieces 
of the (2 3 P |r| 3 S'i) overlap integral. Ebert, Faustov, and 
Galkin [Q include an additional relativistic correction 
to transitions involving the mixed states caused by the 
difference of the c and b quark masses. This leads to 
further differences with the other models. 

Radiative transitions which flip spin are described by 
magnetic dipole (Ml) transitions. The rates for magnetic 
dipole transitions in quarkonium bound states are given 
in the nonrelativistic approximation by [43 . l50| 



r(* - fl) = £ M 2 ^ 3 (2J/ + 1) I (/| Jo(kr/2) I i) | 2 (15) 



TABLE II: El transition rates. The matrix elements were obtained using the wavelunctions ol the GI model 6]. For mixed 
states such as 1P{ and lPi the widths are calculated using the wavefunction conventions defined in eqn. Illl with the mixing 
angles given in Table I and the matrix elements and Cji factors corresponding to the labelling of column 6. 



Initial 


Final 


M % 


M f 


U) 


*-/ 




(f\r\i) 


Cfi 


Width 


state 


state 


(MeV) 


(MeV) 


(MeV) 






(GeV" 1 ) 




(keV) 


1 3 P 2 


l 3 Si 7 


6768 


6338 


416 


1 3 P 2 -> 


l 3 ^! 


1.421 


1 
3 


83 


1P{ 


l 3 Si 7 


6750 


6338 


399 


l 3 Pi -» 


l 3 ^! 


1.435 


1 
3 


11 




l 1 5o 7 




6271 


462 


l 1 ^- 


l'So 


1.288 


1 
3 


80 


lPi 


l 3 S l7 


6741 


6338 


391 


l 3 P 2 -> 


l 3 ^! 


1.435 


1 
3 


60 




l 1 So7 




6271 


454 


l 1 ^- 


l'So 


1.288 


1 
3 


13 


l 3 Po 


l 3 5i 7 


6706 


6338 


358 


1 3 P, -> 


l 3 5i 


1.443 


1 
3 


55 


2 3 Si 


1 3 P 2 7 


6887 


6768 


118 


2 3 5i -> 


1 3 P 2 


-1.914 


5 
9 


5.7 




lPi'7 




6750 


136 


2 3 5i -> 


1 3 P! 


-1.777 


1 
3 


0.7 




IP17 




6741 


144 


2 3 5i -* 


l 3 Pl 


-1.777 


1 
3 


4.7 




l 3 Po7 




6706 


179 


2 3 5'i -> 


l 3 Po 


-1.620 


1 
9 


2.9 


2 1 5'o 


1Pi7 


6855 


6750 


104 


2 1 S -* 


l X Pl 


-1.909 


1 


6.1 




IP17 




6741 


113 


2 1 SW 


l X Pl 


-1.909 


1 


1.3 


1 3 D 3 


1 3 P 2 7 


7045 


6768 


272 


1 3 D 3 ^ 


1 3 P 2 


2.383 


2 


78 


ID' 2 


1 3 P 2 7 


7036 


6768 


263 


1 3 D 2 -H 


• 1 3 P 2 


2.389 


1 
10 


8.8 




lPi'7 




6750 


280 


l 1 ^ -1 


1^1 


2.306 


2 


63 




IP17 




6741 


289 


1 X D 2 -. 


1^1 


2.306 


2 


7.0 


1D 2 


1 3 P 2 7 


7041 


6768 


268 


1 3 D 2 -H 


1 3 P 2 


2.389 


1 
10 


9.6 




lPi'7 




6750 


285 


1 3 D 2 -> 


■ l 3 Pl 


2.274 


3 
10 


15 




IP17 




6741 


294 


1 3 D 2 -. 


• l 3 Pl 


2.274 


3 
10 


64 


l 3 L>i 


1 3 P 2 7 


7028 


6768 


255 


l 3 Di -. 


• 1 3 P 2 


2.391 


1 


1.8 




lPi'7 




6750 


273 


l 3 Di -i 


• l 3 Pl 


2.281 


1 

(i 


4.4 




IP17 




6741 


281 


l 3 Di -i 


• l 3 Pl 


2.281 


1 

(i 


28 




l 3 Po7 




6706 


315 


l 3 Di -i 


■l 3 Po 


2.152 


2 
9 


55 


2 3 P 2 


2 3 Si 7 


7164 


6887 


272 


2 3 P 2 -> 


2 3 Si 


2.195 


1 
3 


55 




l 3 Si 7 




6338 


778 


2 3 P 2 -> 


I 3 Si 


0.2308 


1 
3 


14 




1 3 D3 7 




7045 


118 


2 3 P 2 -> 


1 3 D 3 


-2.072 


14 
20 


6.8 




1^ 2 7 




7036 


127 


2 3 P 2 -> 


l 3 D 2 


-1.970 


1 
10 


0.7 




1L> 2 7 




7041 


122 


2 3 P 2 -> 


1 3 D 2 


-1.970 


1 

10 


0.6 




l 3 £>i7 




7028 


135 


2 3 P 2 -► 


l 3 L>i 


-1.866 


1 

150 


0.1 


2P{ 


2 3 Si 7 


7150 


6887 


258 


2 3 Pi -f 


2 3 S'i 


2.319 


1 

3 


5.5 




2'SoT 




6855 


289 


2 1 P 1 -* 


2 1 5'o 


2.046 


1 
3 


52 




l 3 Si 7 




6338 


769 


2 3 Pi -> 


l 3 Sl 


0.155 


1 
3 


0.6 




1 So 7 




6271 


825 


2 1 Pi -> 


l^So 


0.254 


1 
3 


19 




LD 27 




7036 


113 


2 3 P! -» 


l 3 D 2 


-2.096 


1 
2 


5.5 




1D 2 7 




7041 


108 


2 x Pi -» 


1 1 D 2 


-2.080 


2 
3 


1.3 




l 3 £>i7 




7028 


121 


2 3 Pi ->• 


l 3 L>i 


-1.996 


1 

(i 


0.2 


2Pi 


2 3 Si 7 


7145 


6887 


253 


2 3 P! -» 


2 3 5'i 


2.319 


1 
3 


45 




2'SoT 




6855 


284 


2 x Pi -> 


2 1 S 


2.046 


1 
3 


5.7 




l 3 Si 7 




6338 


761 


2 3 Pi -> 


1 3 S 1 


0.155 


1 
3 


5.4 




1^07 




6271 


820 


2 x Pi -»■ 


I 1 So 


0.254 


1 
3 


2.1 




1^ 2 7 




7036 


108 


2 3 Pi -> 


l 3 D 2 


-2.096 


1 
2 


0.8 




1L> 2 7 




7041 


103 


2 x Pi -> 


l x D 2 


-2.080 


2 
3 


3.6 




l 3 Di7 




7028 


116 


2 3 Pi -> 


l 3 Di 


-1.996 


1 
(i 


1.6 


2 3 Po 


2 3 S l7 


7122 


6887 


231 


2 3 Po^ 


2 3 5' 1 


2.437 


1 
3 


42 




l 3 Si 7 




6338 


741 


2 3 Po^ 


I 3 Si 


0.066 


1 
3 


1.0 




1 3 £>i7 




7028 


93 


2 3 Po -» 


l 3 L>i 


-2.128 


2 
3 


4.2 



> 






7600 



7200 - 



6800 - 



6400 - 



6000 



J572 


7588 
7272 

6887 

6338 
1" 


7455 
7122 

6706 


7475 


7565 7571 


7568- 


.7250 


mi 7W , 


7271" 


'7472 

7150 
'7145 

6750 
"6741 


7365 7376 7380 

7372 7269 7276 


.6855 


7164 7266 
7028 7041 7045 


- 


7036 
6768 

B c Mass Spectrum 


-6271 


- 



lc 3 o 3p p 3 p V) F) V) 3 p p 3 p 

^0 ^1 r r l r 2 U l ^2 U 3 r 2 r 3 V A 



FIG. 1: The B c mass spectrum. 



TABLE III: El transition rates (continued). 



Initial 


Final 


M 8 


M f 


Ul 


»-/ 




(f\r\i) 


C f i 


Width 


state 


state 


(MeV) 


(MeV) 


(MeV) 






(GeV" 1 ) 




(keV) 


1 3 F 4 


1 3 F> 37 


7271 


7045 


222 


1 3 F 4 -» 


1 3 F» 3 


3.156 


3 

7 


81 


1*3 


1 3 F> 37 


7266 


7045 


218 


1 3 F 3 ^ 


1 3 F» 3 


3.159 


1 

21 


3.7 




1F> 27 




7036 


226 


l 1 ^- 


l x D 2 


3.104 


3 

7 


78 




1D 2 7 




7041 


222 


1 3 F 3 ^ 


1 3 D 2 


3.091 


8 

21 


0.5 


1F 3 


1 3 £>3 7 


7276 


7045 


227 


1 3 F 3 -> 


1 3 D 3 


3.159 


1 

21 


5.4 




lD' 2l 




7036 


236 


l 1 ^-* 


l x D 2 


3.104 


3 

7 


0.04 




1D 2 7 




7041 


231 


1 3 F 3 -+ 


1 3 D 2 


3.091 


8 
21 


82 


1 3 F 2 


1 3 F> 37 


7269 


7045 


221 


1 3 F 2 -> 


1 3 D 3 


3.160 


1 

->2o 


0.4 




1£>2 7 




7041 


224 


1 3 F 2 -» 


l 3 D 2 


3.095 


1 
IS 


6.3 




1D 2 7 




7036 


229 


1 3 F 2 -> 


1 3 D 2 


3.095 


1 
15 


6.5 




l 3 F>i 7 




7028 


237 


1 3 F 2 -> 


1 3 D 1 


3.026 


9 

25 


75 



where 



/' 



e c 



'»/-, 



(16) 



e c and e^ are the c-quark and 6-antiquark charges in units 
of |e| (e c = 2/3 and e^ = 1/3), and m c and mj, are the 
quark masses given above. 

The Ml widths and overlap integrals are given in Ta- 
ble El They are compared to other calculations in Table 
IVII Transitions in which the principle quantum number 
changes are refered to as hindered transitions which are 
not allowed in the non-relativistic limit due to the orthog- 



onality of the wavefunctions. Ml transitions, especially 
hindered transitions, are notorious for their sensitivity to 
relativistic corrections [5l| ■ In our calculations the wave- 
function orthogonality is broken by including a smeared 
hyperfine interaction directly in the Hamiltonian so that 
the 3 Si and ^o states have slightly different wavefunc- 
tions. Ebert et al. are more rigorous in how they include 
relativistic corrections [l2j but to improve the J/tp — > rj c ^ 
result they modify the confining potential by making it a 
linear combination of Lorentz vector and Lorentz scalar 
pieces. 



TABLE IV: Comparison of predictions for El transition rates. The column labelled GI is the present work. We quote the 
predicted rates from the various references and do not attempt to normalize the rates to common phase space factors. 



Initial 


Final 






Widths 


(keV) 






state 


state 


GI [6] 


EFG [12] 


GKLT[5] 


EQ [2] 


GJ [9] 


FU [10] 


1 3 P 2 


l 3 Si 7 


83 


107 


102.9 


112.6 


73.6 


126 


1P{ 


l 3 Si 7 


11 


13.6 


8.1 


0.1 


10.5 


26.2 




l'Sbi 


80 


132 


131.1 


56.4 


66.6 


128 


lPi 


l 3 Si 7 


60 


78.9 


77.8 


99.5 


49.0 


75.8 




l 1 So7 


13 


18.4 


11.6 


0.0 


16.6 


32.5 


l 3 Po 


l 3 Si 7 


55 


67.2 


65.3 


79.2 


43.0 


74.2 


2 3 Si 


1 3 P 2 7 


5.7 


5.18 


14.8 


17.7 


4.0 


14.5 




IP17 


0.7 


0.63 


1.0 


0.0 


0.6 


2.5 




IP17 


4.7 


5.05 


12.8 


14.5 


3.6 


13.3 




l 3 Po7 


2.9 


3.78 


7.7 


7.8 


2.6 


9.6 


2 1 So 


1Pi7 


6.1 


3.72 


15.9 


5.2 


3.6 


13.1 




IP17 


1.3 


1.02 


1.9 


0.0 


1.3 


6.4 


1 3 £> 3 


1 3 P 2 7 


78 


102 


76.9 


98.7 






ID* 


1 3 P 2 7 


8.8 


12.8 


6.8 










1Pi7 


63 


116 


46.0 


92.5 a 








IP17 


7.0 


7.25 


25.0 








1D 2 


1 3 P 2 7 


9.6 


27.5 


12.2 


24.7 a 








1Pi7 


15 


14.1 


18.4 


0.1 " 








IP17 


64 


112 


44.6 


88.8 a 






l 3 Di 


1 3 P 2 7 


1.8 


5.52 


2.2 


2.7 








1Pi7 


4.4 


7.66 


3.3 


0.0 








IP17 


28 


73.8 


39.2 


49.3 








l 3 Po7 


55 


128 


79.7 


88.6 






2 3 P 2 


2 3 Si 7 


55 


57.3 


49.4 


73.8 








l 3 ^ 


14 




19.2 


25.8 








l 3 As7 


6.8 


1.59 


10.9 


17.8 








lD' 2l 


0.7 


0.113 


0.5 










I-D2 7 


0.6 


0.269 


1.5 


3.2 a 








l 3 £>i7 


0.1 


0.035 


0.1 


0.2 






2P{ 


2 3 Si 7 


5.5 


9.1 


5.9 


5.4 








2 1 5'o7 


52 


72.5 


58.0 










1 3 S X 7 


0.6 




2.5 


2.1 








1'Sot 


19 




20.1 










1D 27 


5.5 


1.2 


3.5 










I-D2 7 


1.3 


0.149 


2.5 


11.5 a 








l 3 Di7 


0.2 


0.073 


0.3 


0.4 






2Pi 


2 3 Si 7 


45 


37.9 


32.1 


54.3 








2 1 5'o7 


5.7 


11.7 


8.1 










l 3 ^ 


5.4 




15.3 


22.1 








l'5o7 


2.1 




3.1 










1-D 2 7 


0.8 


0.021 


1.2 










I-D27 


3.6 


0.418 


3.9 


9.8 a 








1 3 £>i7 


1.6 


0.184 


1.6 


0.3 






2 3 P> 


2 3 5i7 


42 


29.2 


25.5 


41.2 








l 3 ^^ 


1.0 




16.1 


21.9 








l 3 Di7 


4.2 


0.036 


3.2 


6.9 







a We identify the 1-Dj an d 1-^2 states with the J D 2 and 3 D 2 states 
of Eichten and Quigg 0] . 
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Given the sensitivity of radiative transitions to details 
of the models, precise measurements of electromagnetic 
transition rates would provide stringent tests of the var- 
ious calculations. 



IV. HADRONIC TRANSITIONS 

Hadronic transitions between quarkonium levels are 
needed to estimate branching ratios and discuss search 
strategies for these states. In fact, these are the dominant 
decays for both the ip(2S) and T(25) states. There have 
been numerous theoretical estimates of hadronic transi- 
tions over the years [13, E3, H3, H3, E3, HH ES 113, ES 
153. l60L l6ll l62j |. In some cases the estimates disagree 
by orders of magnitude |54j . Hadronic transitions are 
typically described as a two-step process in which the 
gluons are first emitted from the heavy quarks and then 
recombine into light quarks. A multipole expansion of 
the colour gauge field is employed to describe the emis- 
sion process where the intermediate colour octet quarko- 
nium state is modeled by some sort of quarkonium hybrid 
wavefunction. However, the main disagreement between 
predictions arises from how the rehadronization step is 
estimated. To some extent this latter uncertainty can 
be reduced by employing the multipole expansion of the 
colour g aug e fields developed by Yan and collaborators 
|H IM Hi El together with the Wigner-Eckart theo- 
rem to estimate the El-El transition rates [32| and fix- 
ing the reduced matrix elements by rescaling measured 
transition rates. When no measured transitions exist we 
instead rescale the theoretical estimates of related ma- 
trix elements J33. This is the approach used by Eichten 
and Quigg Q- In addition to El-El transitions there will 
be other transitions such as 3 5i -^ 3 S\ + rj which goes 
via Ml-Ml & E1-M2 multipoles and spin-flip transitions 
such as 3 5i -*- 1 Pi7T7r which goes via El-Ml |2Jj. These 
transitions are suppressed by inverse powers of the quark 
masses and are expected to be small compared to the El- 
El and electromagnetic transitions. The 2 3 S*i — > l 3 S'i+77 
transitions are further suppressed due to being at the 
limit of available phase space. 

The differential rate for El-El transitions from an ini- 
tial quarkonium state $' to the final quarkonium state 
$, and a system of light hadrons, h, is given by the ex- 
pression [32ll33j: 



dT 
~dM< 



[$'^<D + ^(2J+l)]Tj a j j,\ A ^'A) 

(17) 
where £' , £ are the orbital angular momentum and J', J 
are the total angular momentum of the initial and final 
states respectively, s is the spin of the QQ pair, M 2 is the 
invariant mass squared of the light hadron system, { '"} 
is a 6 — j symbol, and Ak(£' , £) are the reduced matrix 
elements. The magnitudes of the Ak(£', £) are model 
dependent with a large variation in their estimates. In 
the soft-pion limit the A\ contributions are suppressed 



so, as is the usual practice, we will take A\{1' ', £) = 0. 
For the remaining rates we use scaling arguments taking 
measured rates as input or, when no measured rates exist, 
we rescale the rates predicted for the bb system by Kuang 
and Yan [33] to obtain the cb rates. The amplitudes for 
El-El transitions depend quadratically on the interquark 
separation so the scaling law between a cb rate and the 
corresponding QQ rate is given by J33 



T(cb) 



(r 2 (cb)} 2 



r(QQ) {r 2 {QQ)Y 



(18) 



up to phase space. The scaling factors used to relate the 
input rates to the cb rates are given in Table IVrTl 

There is a large variation in the reduced rates. For 
example, for the transition l 3 Di — > l 3 5i + tttt estimates 
for ^2(2,0) differ by almost three orders of magnitude 
|33L 154 l6ll I62J] . We point this out as a cautionary 
note to the reader. The reduced rates are summarized 
in Table Eiil For the 2 3 S*i -> l 3 Si + tttt reduced rate 
we take an average of the results taken from rescaling 
the tp' -* J ftp + tttt and T(25) -> T + tttt rates. The 
3 3 S! -> 2 3 S'i +7T7T and 3 3 5i -> 2 3 S 1 +tttt values were ob- 
tained by rescaling the corresponding bb transitions [63J ] . 
A (l, 1) and A 2 (l, 1), corresponding to 2 3 P — > 1 3 Po+tttt 
and 2 3 i-2 — * 1 3 -Pi + tttt respectively, were obtained by 
rescaling the appropriate bb rate predictions given by 
Kuang and Yan J33]. As pointed out above there is con- 
siderable variation in the ^(2,0) amplitude needed for 
the 1 3 -Dj — > l 3 5i +7T7T transitions. The largest predicted 
rate for the Y(1D) — » T(15) + 7T7T transition comes from 
Kuang and Yan l33l and has been ruled out by a re- 
cent CLEO limit Q. The CLEO limit is about a fac- 
tor of 3 larger than the rate predicted by Ko 62] . The 
reduced rate for the cb system found by rescaling the re- 
cent BES measurement 65] of the ip" — > J/tfj-rnr rate is 
considerably larger than the rate found by rescaling the 
T(LD) -> T(15) + tttt CLEO limit. However, it is likely 
that one can reconcile the bb and cc results by properly 
taking into account 2 3 S'i - l 3 Z?i mixing |&1 |g& |£2| . We 
will therefore assume a reduced rate of ^(2, 0) = 21 keV 
which is based on the CLEO limit on the transition 
T(LD) -> T(15) + tttt 64]. The ID -+ IStttt transi- 
tions is the subject of recent interest [5J, |6y, |67| and as 
the experimental measurements improve it would be use- 
ful to revisit these calculations. The uncertainty in these 
hadronic transitions could easily lead to factors of two 
errors in the resulting branching ratios. A final note is 
that we have not considered coupled channel effects to 
DD and BB for the cc and bb states respectively which 
could make a considerable contribution to states close to 
threshold [67|. 

The reduced rates of Table IVIII are used to obtain the 
cb hadronic transitions which are summarized in Table 
IVIIII We do not include decays of the type 2 3 > 1 Pj — > 
1 3 ' 1 P 7 /, as they are expected to be small compared to 
the decays we included. Likewise, transitions with 77 and 
7T° in the final state are possible but are expected to 
have much smaller partial widths and 77 transitions are 



TABLE V: Ml transition rates. The matrix elements were obtained using the wavefunctions ol the GI model 



Initial 


Final 


Mi 


M f 


U) 


(f\j (kr/2)\i) 


Width 


state 


state 


(MeV) 


(MeV) 


(MeV) 


(GeV- 1 ) 


(keV) 


l 3 5i 


l l So 7 


6338 


6271 


67 


0.995 


0.08 


2 3 S' 1 


2 1 5' 7 


6887 


6855 


32 


0.992 


0.01 




l 1 5o 7 




6271 


588 


0.102 


0.6 


2 1 S'o 


l 3 ^ 


6855 


6338 


498 


-0.57 


0.3 


3^! 


^So 7 


7272 


7250 


22 


0.992 


0.003 




2 1 5'o7 




6855 


405 


0.109 


0.2 




l 1 5'o7 




6271 


932 


0.05 


0.6 


^So 


2 3 S l7 


7250 


6887 


354 


-0.04 


0.06 




l 3 Si7 




6338 


855 


0.09 


4.2 



TABLE VI: Comparison ol Ml transition rates. Note that no 
effort has been made to scale the results to a common phase 
space. 



Initial 


Final 


GI 


EFG [12] 


GKLT 


_ 


EQ[2j 


FU [10] 


state 


state 






Widths 


(eV) 




l 3 5i 


l 1 5'o7 


80 


33 


60 




134.5 


59 


2 3 5i 


2 1 So7 


10 


17 


10 




28.9 


12 




l 1 S'o7 


600 


428 


98 




123.4 


122 


2 1 S 


l 3 Si7 


300 


488 


96 




93.3 


139 



further suppressed by phase space. Although the 3 3 Si 
and 3 1 So states are expected to be above BD threshold, 
and therefore relatively broad, we include the two pion 
transitions for the sake of completeness. 



V. WEAK DECAYS 

The final ingredient needed in a study of B c phe- 
nomenology is the B c width and its weak decay partial 
widths. The details of B^ decay have been given else- 
where [H El El El El El El El El E| For com- 
pleteness we give a brief overview of the essential features 
of these decays and summarize the weak decay branch- 
ing ratios in Table E] We refer the interested reader 
to the orig inal literature for details of the calculations 

EH Ea El El El El El El El Ei 

For a rough estimate of the B c width we can treat 
the 6-quark and c-quark decay independently so that 
B c decay can be divided into three classes: (i) the b- 
quark decay with spectator c-quark, (ii) the c-quark de- 
cay with spectator 6-quark, and (iii) the annihilation 
B+ — > i + V£ (cs, us), where £ — e, /i, r. The total 
width is the sum over partial widths 

T{B C -> X) = T(b -> X) + r(c -► X) + T(ann) (19) 

In addition there is a Pauli interference contribution to 
the b — » ccs decay from the c-quark spectator which we 



ignore in this crude estimate. 
In the spectator approximation: 



r(6->A) 



9G 2 F \V cb \ 2 ml 

1927T 3 



and 



r( c ^A) = 



'F\ *C8\ 



\2 rrt b 



1927T 3 



4.8 x 1CT 4 eV 



3.3 x 10" 4 eV 



(20) 



(21) 



where we used \V cb \ = 0.0412, \V CS \ = 0.974, m b = 
4.25 GeV, and m c = 1.25 GeV |p. 

Annihilation widths such as cb — > £v£ are given by the 
expression 



G 



r = if|H c | 2 /iM Bc ^mMi 



M 2 Bc 



a (22) 



where rrii is the mass of the heavier fermion in the given 
decay channel. For lepton channels Cj = 1 while for 
quark channels Cj = 3 1 V qq i \ 2 . The pseudoscalar decay 
constant, /b c , is defined by: 



(Ol&^h^s c(a;)|B c (fc)} = tf Bc V cb k^ 



(23) 



where V cb is the cb element of the Cabibbo-Kobayashi- 
Maskawa matrix, and fc^ is the four-momentum of the 
B c meson. In the non-relativistic limit the pseudoscalar 
decay constant is proportional to the wavefunction at the 
origin and is given by the van Royen-Weisskopf formula 



fB. = ^(0). 



(24) 



This result is modified by QCD corrections and relativis- 
ts effects, which are included using the Mock-Meson ap- 
proach or other relativistic quark models. The predic- 
tions of the various calculations including Lattice QCD 
are summarized in Table llXl Using the mock- meson re- 
sult of /b c = 410 MeV [24|, which is consistent with most 
other predictions, leads to the annihilation width of 



F(onn) = 67 x 10~ 6 eV 



(25) 
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TABLE VII: Estimates of reduced rates for El-El hadronic transitions between cb levels. 



Transition 



(QQ) : rate (keV) 



(r 2 (c6)>/{r 2 (QQ)> 



Reduced cb rate (keV) 



2 3 Sl -> l 3 5l + 7T7T 


(cc) 


146 ± 14 a 




(66): 12.2 ±2 " 




Average 


3 3 Sl -> 2 3 Sl + 7T7T 


(66) 


1.26 ±0.25 a 


3 3 5l -» l 3 5l + 7T7T 


(66) 


1.72 ±0.25 a 


2 3 P -> 1 3 P() + 7T7T 


(66) 


0.4 b 


2 3 P 2 _» 1 3 P! + 7T7T 


(66) 


0.01 b 


l 3 Dl -> l 3 Sl +7T7T 


(cc) 


120 ± 57 c 




(cc) 


< 92 90% C.L 




(66) 


24 e 




(66) 


0.07 / 




(66) 


0.56 ± 0.07 9 



"From PDG Ref. El. 

'From Kuang and Yan using model C: modified Richardson po- 
tential Ref. l33l . 

c From BES R.ef.|6fij. 

d CLEO 90% C.L. upper limit Id^ . 

c From Kuang and Yan using model A: linear plus Coulomb po- 
tential Ref. J33J . 

^From Moxhay Ref. |6lj . 

9 From BR prediction of Ko l62l combined with the total width 
predicted by Kuang and Rosncr 45]. 



0.75 
1.63 

1.56 
1.56 
1.56 
1.57 
0.78 
0.78 
1.6 
1.6 
1.6 



A o (0,0) = 82 ±8 

Ao(0,0) = 33 ±5 

57 ±7 

A' o (0,0) =3.1 ±0.6 

Aq(0,0) =4.2 ±0.6 

Ao(l,l) = 2.92 

A 2 (l,l) = 0.164 

A 2 (2,0) = 360 ±170 

A 2 (2,0) <280 

A 2 (2,0) = 307 

A 2 (2,0) =0.9 

A 2 {2,0) = 7.2 



Adding this result to the spectator contributions gives 
T{total) — 8.8 x 10~ 4 eV corresponding to a B c life- 
time of t = 0.75 ps which is in rough agreement with 
the measured value of r = 0.46l^i6 P s - A more care- 
ful calculation by Kiselev gives r = 0.5 ps 69]. The 
approximate branching fractions for the 6-decay, c-decay, 
and annihilation processes are 54%, 38%, and 8% respec- 
tively. These BR's are modified by strong interaction ef- 
fects H H |zl III III IH which are included in some 
recent calculations of BR's. The branching ratios of some 
prominent decay modes are summarized in Table IXl In 
addition, weak B c decays to P-wave charmonium states, 
Xc or h c , are potentially important decay modes |76j but 
we will neglect them in favour of simpler to observe decay 
chains. 



VI. EXPERIMENTAL SIGNATURES AND 
SEARCH STRATEGIES 

B c mesons offer a rich spectroscopy of narrow states 
to study; there are 2 S'-wave, 2 P-wave, and 1 D-wave 
B c multiplets below BD threshold. Because B c mesons 
carry flavour they cannot annihilate into gluons and are 
expected to be quite narrow; < 100 keV. In addition, the 
F-wave states are just above threshold so might also be 
relatively narrow due to the angular momentum barrier 
which would suppress the decay [79j. Two ingredients 
are necessary for the study of B c spectroscopy; that they 
be produced in sufficient quantity and that they yield a 
signal that can be distinguished from background. 



B c production production proceeds via the hard as- 
sociative production of the two heavy quark pairs cc 
and bb which suppresses the B c yield relative to beauty 
hadrons by O(10 -3 ) [21|. Because fragmentation dom- 
inates at high-pT it has proven useful to describe B c 
meson production by hadronization of individual high- 
Pt partons using the factorization formalism based on 
non-relativistic QCD [13. This approach was utilized 
by a number of authors to calculate B c production at 
hadron colliders. At the Tevatron, with acceptance cuts 
of p t > 6 GeV and \y\ < 1, Cheung estimates that O(10 7 ) 
B c mesons should be produced for an integrated lumi- 
nosity of 1 fb" 1 [23. At the LHC O(10 9 ) B c mesons 
are expected to be produced for 100 fb -1 with the kine- 
matic cuts oi p t > 10 GeV and \y(B c )\ < 2.5 22]. The 
B* cross sections are expected to be 50-100% larger than 
the B c cross sections [l^, |2l|, |22| . The pr distribution 
falls rapidly so that the small px region is quite impor- 
tant [l|| but also making the cross sections sensitive to 
the exact values of the kinematic cuts [nj, |2(j ■ With the 
high luminosity of the LHC one expects a sizable num- 
ber of P- and D-wave (cb) states (~ 20% and ~ 2% of 
the total inclusive B c cross section respectively) as well 
as excited 5-wave states to be produced (2S'/1S' ~ 0.6) 
miHEl. The LHC should therefore produce sufficient 
B c mesons to allow the study of the cb spectroscopy and 
decays. We will use these numbers as the starting point 
to estimate the number of J5 c 's produced in a particular 
decay chain. 

Irregardless of which state is produced, it will even- 
tually cascade decay to the B c ground state via electro- 
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TABLE VIII: Rates for two-pion El-El hadronic transitions. The reduced rates are denoted by Ak(i' ', tj where k is the rank 
of the irreducible tensor and £' and £ are the orbital anular momenta of the initial and final states. 

Transition Expression for Rate (c6) rate (keV) 

A (0, 0) 57 ± 7 

A o {0,0) 57 ±7 

A' {0, 0) 3.1 ±0.6 

A' {0, 0) 3.1 ±0.6 

A'o(0, 0) 4.2 ± 0.6 

Aq(0, 0) 4.2 ±0.6 

2 3 P 2 -* 1 3 P 2 + tttt |A (1, l) + i J 4i(l, 1) + |A 2 (1, 1) 1.0 

2 3 P 2 -> IP/ ± 7T7T ±Ai(l, 1) + -kA a (l, 1) a 0.004 6 

2 3 P 2 -> lPi + 7T7T iil(l,l) + |4 2 (l,l)° 0.021 6 



2 3 5i - 


-> l 3 ^! ± 7T7T 


2 1 S , - 


-> l 1 So + tttt 


3 3 S! - 


-> 2 3 5l ± 7T7T 


S^o - 


■>2 1 S + irir 


S 3 ^! - 


->• l 3 Sl ± 7T7T 


3^0- 


-> 1 So+ TTTT 



2 3 P 2 -f 1 3 P ± TTTT T5^2(l, 1) 0.011 

2P{ -» 1 3 P 2 ± 7T7T |Al(l, 1) + ±A 2 (1, 1) c 0.004 ' 

2Pi -> 1 3 P 2 ± tttt |4i(l, 1) + jA 2 (l, 1) c 0.037 b 

2P[^IP[+tttt A (l, 1) + Ai(l, l) + iA 2 (l, 1) d 1.2 e 

2PJ 1 -> lPi ± tttt A (l, l) + vli(l, l) + |A a (l, 1) d 0.1 c 

2Pl -> IP; + TTTT |A (1, l)+i^i(l, l) + i J 4 2 (l, 1)' 0.02 e 

2Pi^lP 1 ±7T7r |A (1, l) + ^Ai(l, 1) + ±A 2 (1, l) s 2.7 e 

2P^ -» 1 3 P ± 7T7T |Al (1, 1) 9 

2Pl -> 1 3 P ± 7T7T |Al(l, l) 9 

2 3 P -> 1 3 P 2 ± tttt |-4a(l, 1) 0.0547 

2 3 P -> lPi ± 7T7T |Al(l, 1) ' l 

2 3 P -> lPl ± 7T7T |^4l(l, 1) A 

2 3 P -> 1 3 P + 7T7T ly4 (l, 1) 0.97 

l 3 _Dl,3 -» l 3 Si ± tttt M2(2, 0) * 4.3 

1D' 2 -> 1 3 S*1 ± TTTT 1,42(2, 0) ! 2.1 6 

1D 2 -» l 3 Si ± 7T7T |A 2 (2, 0) ' 2.2 b 

1D' 2 -» 1^0 ± 7T7T 1-42(2, 0) * 2.2 b 

\D 2 -> l 1 So±7r 7 r 1^2(2,0)* 2.1 6 

"The expression is for the 3 P 2 — > 3 Pi transition. 

''These rates include the appropriate mixing angles defined in eqn. 
and given in Table |l] 

The expression is for the 3 Pi — > 3 P 2 transition. 

d The expression is for the X P\ -^ 1 P± transition. 

These rates include the appropriate mixing angles. We assume 
a +ve phase between the 1 Pi — » J Pi and 3 Pi — > 3 Pi amplitudes. 
The 2Pj — > lPi and 2Pi — * 1P^ widths are most sensitive to this 
phase. In any case the widths are expected to be quite small and 
not particularly important. 

The expression is for the 3 Pi — » 3 Pi transition. 

The expression is for the 3 Pi — > 3 Pq transition. 

The expression is for the 3 Po — » 3 Pi transition. 

'These correspond to the 3 Dj — > 3 Si or 1 D 2 — ► 1 So transitions 
as appropriate. 





TABLE IX: Comparison of predictions for the pseudoscalar decay constant of the B c meson. 




GI f34] 


EFG [12] GKLT [5] EQ [2] Fu [10] Ki[70] Lattice [44] 


Lattice [711 


410 ± 40 


433 500 ± 80 500 a 517 395 ± 15 440 ± 20 


420 ± 13 



"Using Buchmiiller-Tye potential. 
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TABLE X: Branching ratios in % for some prominent exclu- 
sive B+ decays. 



Mode 


Kis [69] 


IKS f72] 


CC ° 


f73] EFG [75] 


NW [74] 


Bt - 


J/ipe + v 


1.9 


2.3 


2.4 


1.2 


1.5 


-> 


r) c e + v 


0.75 


0.98 


1.0 


0.42 


0.52 


-> 


B°e+v 


0.34 


0.15 


0.16 


0.042 


0.05 


-» 


B*°e+v 


0.58 


0.16 


0.23 


0.12 


0.05 


-> 


B° s e+u 


4.03 


2.0 


1.9 


0.84 


0.94 


-> 


Bfe+u 


5.06 


2.6 


3.1 


1.75 


1.44 


-» 


J/t/jTT + 


0.13 




0.22 


0.061 




-> 


J/^ P + 


0.40 




0.66 


0.16 




-» 


77 c 7T + 


0.20 




0.23 


0.085 




-> 


?7cP + 


0.42 




0.61 


0.21 




-> 


B°7r+ 


16.4 




5.1 


2.52 




-» 


B.V 


7.2 




3.9 


1.41 




-> 


B S *V+ 


6.5 




4.5 


1.61 




-> 


Bfp + 


20.2 




13.1 


11.1 




-> 


B° 3 K+ 


1.06 




0.37 


0.21 




-> 


BfK+ 


0.37 




0.26 


0.11 




-> 


B°7V + 


1.06 




0.29 


0.10 




-» 


B°p+ 


0.96 




0.52 


0.13 




-> 


B*°TV + 


0.95 




0.25 


0.03 




-» 


B*°p+ 


2.57 




1.0 


0.68 




-» 


B+K° 


1.98 




0.30 


0.24 




-> 


B+K*° 


0.43 




0.21 


0.09 




-» 


B'+K° 


1.60 




0.23 


0.11 




-» 


B*+K*° 


1.67 




0.44 


0.84 




~ * 


T + V T 

cs 


1.6 
4.9 











"Using the PDG 63J central value for the total width of T 
f .43 X fCT 3 eV to obtain the BR's. 



magnetic and hadronic transitions so that the B c must 
be observed in order to reconstruct the parent particle 
for the particular decay chain. Prominent decays of the 
B c are given in Table E] B c decays with a J/ip in the 
final state such as J/ip + X where X can be a 7r + , p + , 
or l + vi are especially useful as the J/ip — > l + £~ pro- 
vides a useful trigger for B c events. The golden chan- 
nel to detect B c is B c — > J/tpir(p) but their BR's are 
quite small O(0.2 — 0.4%) resulting in a combined BR 
for B c -> J/ip + 7T+ -> I'- (! +7T+ of - 0.02%. This 
would yield about 2000 events for 1 fb _1 integrated lu- 
minosity at the Tevatron and about 2 x 10 5 events for 
100 fb _1 integrated luminosity at the LHC. The decay 
B^ — > J/ipi^Pc with J/ip — > e + e~ has a distinctive sig- 
nature of 3 charged leptons coming from a common sec- 
ondary vertex and has a BR of about 2% resulting in a 
combined BR for B c -> J/ip£+v e -» i'-l'+l+vi (£ = e, fj) 
of about 0.2% yielding 2 x 10 4 events at the Tevatron and 
2 x 10 6 at the LHC. The trade off is that the semilcp- 



tonic decay mode B c — > J/tp + tvt has a larger branching 
ratio but also has missing energy while B^ — > J/ipir 
has a smaller BR but has the advantage that the B c can 
be fully reconstructed. In these B c decays the b quark 
decays to charm. B c decays in which the c quark de- 
cays, such as -B+ — > B ( *ir + , have much larger BR's and 
could also prove to be important modes if B c 's can be 
reconstructed in these channels. 

To estimate event rates we also need to include detec- 
tion efficiencies. Simulations by DO and CDF [3J find ef- 
ficiencies for the exclusive decays B c — > J/ijjtt — » n + [i~ n 
and B c — > J/ipli^i — > fi + prtvt of ~ 2% and — 4% re- 
spectively. 

In Table IXII we combine the electromagnectic transi- 
tions widths with the hadronic transition widths to give 
total widths and BR's. The are used in Table IXIIII to 
give estimates for the number of events expected at the 
Tevatron and LHC for the more prominent decay chains 
of B c excited states. We assume that the ground B c 
state is observed in the B c — ► J/ip£ + ve — ► l'~£' + £ + vi 
and B c -> J/ip + n + -»• £'-£' + ir + decay modes with BR's 
2% and 0.2% respectively and detection efficiencies of 
~ 2% and ~ 4% respectively. We include a factor of 2 to 
take into account both the e + e~ and fi + pT decay modes 
of the J/tjj and a factor of 2 to take into account the 
production of both charge conjugate B c states. To take 
into account the relative production rates of different ex- 
cited states we use the production rates given above of 
10 7 B c 's at the Tevatron and 10 9 S c 's at the LHC with 
relative numbers of x2 for the B*, 0.2 for the IP states, 
0.02 for the ID states, and the 0.6 for the 2S relative to 
the IS states. For the 2P states we use the same factor 
of 0.6 for 2P relative to IP but this is a rather arbitrary 
assumption. As noted already, the cross sections are very 
sensitive to the kinematic cuts so the number of events 
expected should only be taken as rough estimates. 

The signal for excited B c states is a photon or pions 
in coincidence with B c decay. A serious omission from 
the estimates given in Table IXIIII is the neglect of tag- 
ging efficiencies for the photons and pions in the transi- 
tions. The photon ID can be relatively high depending 
on the kinematics. For the pions, the combinatorial back- 
ground is large and tt/K separation is not so good so one 
really needs to do studies of specific processes. Thus, 
understanding photon and pion identification requires a 
detailed simulation study which is beyond the scope of 
this paper. However, one might be optimistic, given the 
success of CDF in studying \c Producti on 1801 and the 
observation of the A"(3872) — > J/tpir + ir~ 81] in pp colli- 
sions at the Tevatron. 

Notwithstanding the previous caveat, it should be pos- 
sible to observe the IS, 25, and IP states at the Teva- 
tron. It is also possible, although only marginally so, 
that some of the ID and 2P might also be seen. With 
the higher statistics available at the LHC, all cb states 
below threshold could potentially be observed, although 
the larger backgrounds will make this quite challenging. 

The fact that the B c is not an eigenstate of charge con- 



13 



jugation helps simplify the search for states such as the 
1P{, 1-Pi, ID' 2 , and 1£>2- The singlet component of these 
states allows El or hadronic transitions directly to the 
ground state B c with large BR's. This should simplify 
the reconstruction efforts significantly. For example 1P{ 
production and decay to the B c with its subsequent de- 
cay to J/tpn and J j^ivi should produce 0(600) events in 
Run II at the Tevatron. Likewise, 2P[ and its subsequent 
decay should produce 0(100) events. These yields would 
be enhanced if other B c decay modes with larger BR's 
could be utilized. The discovery of the 2P[ would yield 
important spectroscopic information in addition to being 
an experimental tour de force. Of course, these BR's are 
highly sensitive to the \P[ — 1P\ mixing angle. With 
enough measurements these details can be constrained 
and different models can be differentiated. There will 
also be large number of events for decay chains going 
via an intermediate B* such as 2 3 S\ — > l 3 5i — > B c and 

1 3 P 2 — » I 3 Si — > B c . However, it will be crucial to detect 
the 67 keV photon in B* — ► B c + 7 in these cases, a very 
challenging experimental task. 



model. For the most part the mass predictions are con- 
sistent with other models, within the accuracy of these 
models. The largest discrepancy in predictions is for the 
triplet-singlet mixing angles. This has implications for 
transitions between states so can be tested with appro- 
priate measurements. Combining the BR's we have cal- 
culated with B c production cross sections from the litera- 
ture we see that the IS, IP, and 25 states should be pro- 
duced in sufficient numbers to be observed at the Teva- 
tron. With the higher statistics of the LHC, it should 
also be possible to observe the ID and 2P states. It will 
be a significant experimental challenge to extract the sig- 
nals for these states from the large background but their 
observation would add considerably to our knowledge of 
quarkonium spectroscopy and discriminate between the 
various models that exist in the literature. 
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TABLE XI: Partial widths and branching fractions for strong 
and electromagnetic transitions. Details of the calculations 
are given in the text. 



TABLE XII: Partial widths and branching fractions for strong 
and electromagnetic transitions (continued). 



Initial 


Final 


Width 


B.F. 


state 


state 


(keV) 


(%) 


l 3 Si 


l 1 ^ 


0.08 


100 


1 3 P 2 


l 3 ^! 


83 


100 


1P{ 


l 3 ^! 


11 


12.1 




l 1 ^ 


80 


87.9 




Total 


91 


100 


lPi 


l 3 ^! 


60 


82.2 




l'So 


13 


17.8 




Total 


73 


100 


l 3 Po 


l 3 ^! 


55 


100 



2 1 So 


l^O+TTTT 


57 ±7 


88.1 




lPi+1 


6.1 


9.4 




lPi+7 


1.3 


2.0 




l 3 Si + 7 


0.3 


0.5 




Total 


64.7 


100 


2 3 Si 


l 3 5"l + 7T7T 


57 ±7 


79.6 




1 3 P 2 + 7 


5.7 


8.0 




LPi+7 


0.7 


1.0 




lPi+7 


4.7 


6.6 




l 3 Po + 7 


2.9 


4.0 




2 1 So + 7 


0.01 


1 x 10~ 2 




l 1 S () + 7 


0.6 


0.8 




Total 


71.6 


100 



1 3 P>3 


l 3 5"l + 7T7T 4.3 


5.2 




1 3 P 2 + 7 78 


94.8 




Total 82.3 


100 



IDa 


l 3 ^! + 7T7T 


2.1 


2.5 






l^O+TTTT 


2.2 


2.6 






1 3 P 2 + 7 


8.8 


10.6 






1P/+7 


63 


75.8 






lPi+7 


7 


8.4 






Total 


83.1 


100 




1D 2 


l 3 ^! + 7T7T 


2.2 


2.4 






^So+TTTV 


2.1 


2.3 






1 3 P 2 + 7 


9.6 


10.3 






1 A' + 7 


15 


16.1 


[30] 




lPi+7 


64 


68.9 


[31] 




Total 


92.9 


100 


[32] 
[33] 


l 3 Di 


l 3 ^! + 7T7T 


4.3 


4.6 




1 3 P 2 + 7 


1.8 


1.9 


[34] 




LPi+7 

lPi+7 


4.4 

28 


4.7 
29.9 


[35] 




l 3 Po + 7 


55 


58.8 


[36] 




Total 


93.5 


100 












[37] 



Initial 


Final 


Width 


B.F. 


state 


state 


(keV) 


(%) 


2 3 P 2 


1 3 P 2 + 7T7T 


1.0 


1.3 




2 3 Si+ 7 


55 


70.3 




l 3 5i+ 7 


14 


17.9 




l 3 £>3+7 


6.8 


8.7 




lZ>2+7 


0.7 


0.9 




1D 2 +j 


0.6 


0.8 




l 3 £>i+7 


0.1 


0.1 




Total 


78.2 


100 


2p; 


1P{ + nn 


1.2 


1.4 




2 3 Si+ 7 


5.5 


6.4 




2 1 5 ( ,+7 


52 


61.0 




l 3 5i+7 


0.6 


0.7 




l 1 5o+7 


19 


22.3 




LD 2 +7 


5.5 


6.4 




lP> 2 +7 


1.3 


1.5 




l 3 Di+7 


0.2 


0.2 




Total 


85.3 


100 


2Pi 


lPl + 7T7T 


2.7 


4.0 




2 3 Sx+ 7 


45 


67.3 




2^0+7 


5.7 


8.5 




l 3 5!+ 7 


5.4 


8.1 




l 1 So+7 


2.1 


3.1 




1^2+7 


0.8 


1.2 




1^2+7 


3.6 


5.4 




l 3 Di+7 


1.6 


2.4 




Total 


66.9 


100 


2 3 P 


1 3 P + 7T7T 


1.0 


2.1 




2 3 Si+ 7 


42 


87.1 




l 3 5i+ 7 


1.0 


2.1 




l 3 £>i+7 


4.2 


8.7 




Total 


48.2 


100 
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TABLE XIII: Expected event rates for various decay chains 
at the Tevatron and the LHC. The B c is assumed to decay to 

B c — > J/ipTT — > fi + n~-K and B c — > J/tjj(.V{: — ♦ fi + n~ £ve final 
states. Details of these estimates are described in the text. 



[39] 
[40] 
[41] 

[42] 

[43] 

[44] 

[45] 

[46] 
[47] 
[48] 
[49] 



Decay Chain 


Tevatron 


LHC 




l 3 5i -2* P c 


3.4 x 10 3 


3.4 x 


10 5 


[51] 


Z Do — > D c 


1.8 x 10 3 


1.8 x 


10 5 




2 3 Si ^ r% -i P c 


1.6 x 10 3 


1.6 x 


10 s 




2 3 5l -^ Be 


16 


1.6 x 


10 3 




1 3 P 2 -^ l 3 5i ^ Be 


6.7 x 10 2 


6.7 x 


10 4 




IPi ^ Be 


5.9 x 10 2 


5.9 x 


10 4 




IPi -^ l 3 ^ -^ B c 


5.5 x 10 2 


5.5 x 


10 4 


[52] 


IPi ^ P c 

1 3 P -^ l 3 5l -^ Be 


1.2 x 10 2 
6.7 x 10 2 


1.2 x 
6.7 x 


10 4 

10 4 


[53] 
[54] 


1 3 D 3 ™ l 3 Sx 1 Be 


3.5 


3.5 x 


10 2 




1 3 D 3 -^ 1 3 P 2 -^ l 3 ^! 


-^ Be 64 


6.4 x 


10 3 


[55] 


1P 2 =5 i^o 


1.7 


1.7 x 


10 2 


[56] 
[57] 
[58] 


1P 2 -^ lP^ -^ P c 


45 
1.4 


4.5 x 
1.4 x 


10 3 

If) 2 


1P 2 -^ IPi -^ P c 


9.5 


9.5 x 


10 2 


[59] 


1P 2 -^ IPi -i P c 


8.2 


8.2 x 


10 2 


[60] 


l 3 Px =3 1 3 S! -i P c 


3.1 


3.1 x 


10 2 


l 3 Pl -^ IPi' -^ P c 


2.8 


2.8 x 


10 2 


[61] 


l 3 Pi -^ IPi -^ P c 


3.6 


3.6 x 


10 2 


[62] 


2 3 P 2 ^> 2 3 S' 1 -^ P c 


2.3 


2.3 x 


If) 2 


[63] 


2 3 P 2 ^> l 3 5i -^ P c 


72 


7.2 x 


10 3 


[64] 


2Pi -^ P c 


90 


9.0 x 


10 3 


2Pi ^ P c 


12 


1.2 x 


10 3 


[65] 


2Pi -i 2 3 5'i -^ P c 


2.2 


2.2 x 


If) 2 


[66] 


2 3 P -^ 2 3 5i -2> P c 


2.8 


2.8 x 


10 2 


[67] 
[681 
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